Angular Momentum Transfer in Vela-like Pulsar Glitches 

Pierre M. Pizzochero 

Dipartimento di Fisica, Universitd degli Studi di Milano, and Istituto Nazionale di Fisica 
Nucleare, sezione di Milano, Via Celoria 16, 20133 Milano, Italy 



The angular momentum transfer associated to Vela-like glitches has never 
been calculated directly within a realistic scenario for the storage and release 
of superfluid vorticity; therefore, the explanation of giant glitches in terms of 
vortices has not yet been tested against observations. We present the first phys- 
ically reasonable model, both at the microscopic and macroscopic level (spher- 
ical geometry, n=l polytropic density profile, density-dependent pinning forces 
compatible with vortex rigidity), to determine where in the star the vorticity is 
pinned, how much of it, and for how long. For standard neutron star parameters 
(M = IAMq.Rs = 10 km, (l = f^vcia = -10"^° Hz s"^), we find that maxi- 
mum pinning forces of order fm ~ 10^^ dyn cm~^ can accumulate ALgi ^ 10^° 
erg s of superfiuid angular momentum, and release it to the crust at intervals 
Atgi ~ 3 years. This estimate of ALgi is one order of magnitude smaller than 
what implied indirectly by current models for post-glitch recovery, where the core 
and inner-crust vortices are taken as physically disconnected; yet, it successfully 
yields the magnitudes observed in recent Vela glitches for both jump parameters, 
Afigi and Afigi, provided one assumes that only a small fraction (< 10%) of the 
total star vorticity is coupled to the crust on the short timescale of a glitch. This 
is reasonable in our approach, where no layer of normal matter exists between the 
core and the inner-crust, as indicated by existing microscopic calculation. The 
new scenario presented here is nonetheless compatible with current post-glitch 
models. 



Subject headings: stars: neutron — dense matter — pulsars: general — pulsars: 
individual (PSR J0835-4510) 



Glitches are sudden spin-ups observed in the otherwise decreasing rotational frequency 
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Introduction 




Their origin is still debated: the giant spin-ups observed in the 
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twenty known Vela-like glitchers (Espinoza et al. 2011) could indicate the presence of bulk 



superfluidity inside these stars. In this scenario, giant glitches would represent the natural 
macroscopic outcome of the interaction between quantized neutron vortex lines, which carry 
the angular momentum of the rotating chargeless superfluid, and the Coulomb lattice of 
neutron-rich nuclear clusters, which coexists with the neutron superfluid in the inner crust 



(Negele & Vautherin 1973 ). Indeed, this interaction can pin vortices to the normal component 



of the star, thus freezing the superfluid vorticity and storing its angular momentum. Only 
when the hydrodynamical lift on a vortex (Magnus force), which increases as the pulsar slows 
down, equals the pinning force on a line, the vortex is unbound from the lattice: thus free to 
move under the action of drag forces, it can transfer its angular momentum to the normal 



component of the star. According to Anderson & Itoh (1975), giant glitches are due to the 



sudden and simultaneous depinning of a large number of accumulated vortices, followed by 
the rapid transfer of their angular momentum to the observable normal crust (which consists 
of the outer crust plus all the other charged components in the star, electrons, protons, and 
nuclear clusters, strongly coupled together by the pulsar magnetic fleld). Such a storage and 



trigger mechanism would have a natural periodicity, as indeed observed in Vela (Dodson et 



al. 2007). 



The vortex scenario for glitches was roughly compared to existing observations in a 
simple but instructive toy-model, which assumed a cylindrical, uniform-density star with 



a cylindrical pinning shell (corresponding to the inner-crust) close to the surface (Pines et 



al. 1980 Alpar et al. 1981; Anderson et al. 1982). Although the naive treatment of the 



vortex-nucleus and vortex-lattice interactions gave pinning forces three orders of magnitude 
larger than what required to explain the average interval between glitches observed in Vela 
(Atgi 3 years), the model predicted the correct orders of magnitude for the typical glitch 
parameters known at the time, namely the jump in angular velocity, Afigi ^ 10~^Q, and the 
jump in angular acceleration, AQgi ^ 10~^fi (the pre-glitch, steady-state parameters for Vela 
being flyeia = 70 Hz and fiyeia = —9.8 x 10^^^ Hz s~^). In spite of these positive preliminary 



results. Pines et al. (1980) carefully pointed out that the effects of some crucial corrections 



had to be taken into account before drawing any conclusion: the spherical geometry of the 
star, the radial density proflle required by gravitational equilibrium, the density dependence 
of the pinning interaction, and the presence of different superfluid phases along the star 
proflle. To date, however, this has not been done in any coherent and consistent model; thus, 
the explanation of giant glitches in terms of vortices is not yet tested against observations, 
leaving the origin of these spin-ups a still open question. 

The post-glitch recovery of pulsars, on the other hand, has been successfully interpreted 



in terms of vortex motion under drag forces. Early on, the phenomenological model of Baym 



et al. ( 1969 ) explained the slow relaxation to steady-state following a glitch as due to the weak 
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interaction between a normal and a superfluid component, each rotating rigidly. Following 
a glitch, the response of the model is linear in the initial perturbation, and relaxes back to 
steady-state conditions exponentially, with a relaxation time which is inversely proportional 
to the strength of the interaction between the two components. The simple two-components 
model was then reformulated in terms of vortex motion, to allow for differential rotation of 



the superfluid (Pines et al. 1980 ). Eventually, two scenarios were developed to describe vortex 



dynamics between glitches: thermally- activated creep of strongly pinned vortices (Alpar et 



al. 


1984a|b 


Link et al. 


1993) 



(Jones 1990 1991, 1993). The vortex creep model was motivated by the large pinning forces 



obtained in early calculations (Alpar 1977 Epstein & Baym 1988). Later on, however, the 



microscopic vortex-nucleus interaction was shown to be one order of magnitude smaller than 



what found earlier (Donati & Pizzochero 2003 2004 2006). Moreover, Jones (1992) argued 



that the mesoscopic vortex-lattice interaction, necessary to calculate the macroscopic pinning 
force on a vortex line, is likely to be a factor ai ~ 10~^ smaller than what naively assumed 
in early calculations, due to the random orientation of the macro-crystals forming the inner 
crust, as well as to the rigidity of vortex lines on distances of order 10^ — lO^-R^s (with R^g the 
radius of the Wigner-Seitz cells describing the nuclear lattice). Significantly smaller pinning 
forces favor the corotation model, where unpinned vortex lines are weakly coupled to the 
normal crust by small drag forces; thence, in steady-state they (nearly) corotate with the 
superfluid, while the response to perturbations is linear. Finally, the Christmas 1988 Vela 



glitch (Flanagan 1990) showed that the creep model can fit observations only in the linear 



response regime (Alpar et al. 1989), in which case it is equivalent to the corotation model, 
but with a different temperature-dependence of the drag parameters. Observationally, the 
post-glitch recovery of Vela is well described by a sum of exponential terms with different 
amplitudes and relaxation times ( Dodson et al.|[2002 ), and this can be explained in terms of 
linear response of regions of the superfluid characterized by different drag parameters. The 
dissipative force has been evaluated for several densities of interest, and the corresponding 
drag parameters yield relaxation times and glitch rise-times compatible with observations 



(Jones 1990, 1992 Epstein & Baym 1992) 



Three aspects of current post-glitch models are relevant here: 



i) although simulations successfully reproduce the observed recovery of Vela (Larson & Link 



2002), the glitch itself is always introduced by hand, as an ad hoc initial condition. 

ii) core and crust vortices are taken as physically disconnected, namely a layer of normal 

matter is assumed between the S-wave neutron superfluid found at subnuclear densities and 



the P-wave one found above nuclear saturation (Jones 1990). Recent microscopic calcula- 
tions, however, do not show any discontinuity of this kind ( Zhou et "aL]|2004 ) , thus indicating 
continuous vortex lines throughout the star. 
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iii) in the core, the superconducting state of protons determines in part the mutual friction 
on neutron vortex hues. A type II superconductor corresponds to the strong-drag limit, with 
vortices entangled in a dense array of magnetic flux tubes (Link 2003); for a type I super- 



conductor, instead, both the weak- and strong-drag limits have been suggested (Sedrakian 



2005||Jones|[2006| ). To date, the microscopic nature of the proton superconductor is far from 



settled; therefore, both scenarios of weak and strong mutual friction in the core should be 
taken into account in the study of glitches. Post-glitch models, however, are not affected 
by this theoretical uncertainty, since they involve mostly vortices in the equatorial regions, 
lying entirely in the inner crust. 



1045 g 



while its 
tot with the 



A typical neutron star has total momentum of inertia /tot 
inner crust has lie ~ 10~^/tot- The accidental coincidence of the ratio /ic//i 
early observations of AQg\ ^ lO^^Q and with fact that about 1.7% of the Vela spin-down is 
reversed during a glitch, led to consider glitches as related to crust vorticity alone, and thence 
to the assumption of disconnected neutron superfiuids. This scenario, however, has direct 
implications on the glitch energetics. Indeed, since vortex lines in the core are strongly 



coupled to the normal component, being magnetized by entrainment effects (Alpar et al. 



1984c), the normal crust comprises most of the star and Jc 



'tot- 



This implies that the 



angular momentum transferred during the glitch is ALgi = IcAflgi ^ 10 erg s, and the 



corresponding glitch energy is Ai^gi = ALgiQ ^ 10 erg; both values appear too large. 
On the one hand, 10^^ erg s corresponds to the difference in angular momentum of the 
entire inner crust between glitches, thus requiring some unlikely mechanism that freezes the 
vorticity everywhere in the crust for about 3 years, and then releases it simultaneously. On 
the other hand, observations of the power wind nebula surrounding Vela indicate an upper 
limit of ~ 10^^ erg to the glitch energy (Helfand et al. 2001). 

In this letter we present the first realistic (with several approximations, but still preserv- 
ing the essential physics) and consistent model to determine where in the star the vorticity 
is pinned, how much of it, and for how long. The model has been tested against obser- 
vations using realistic equations of state (EoS) for dense matter and implementing general 
relativistic hydrostatic equilibrium (Pizzochero, Seveso & Haskell, in preparation). More- 



over, initial dynamical simulations based on the multifiuid formalism of Andersson fc Comer 



(2006) confirm the main assumptions and predictions of the model (Haskell et al. 2011). 
Here, however, we will discuss a fully analytical, Newtonian version of the model: it yields 
the correct orders of magnitude for all relevant variables, but provides deeper insight than 
any numerical treatment. 
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2. The model 

We now outline the main assumptions of the model and present the resulting equations; 
details and calculations, together with a parameter study of the solutions, will be given in a 
longer article (Pizzochero, in preparation; from now on. Paper I). 

We describe the core and inner crust as a n=l polytrope, of mass M and radius Rs. The 
actual radius of the star will be larger, because of the overlying outer crust; its presence, 
however, can be ignored here, since it contributes negligibly to the mass and moment of 
inertia of the normal component. The polytropic relation P oc is a very soft EoS for 
dense matter; realistic soft EoSs yield Rs ~ 10 km for M = IAMq. The density profile 
is M = sin(7r^)/(7r^), with the dimensionless radius ^ = r/Rg and the density u = p/X, 
normalized to its central value A = 7rM/(4i?^). The radius of the core, Rc, corresponds to 
the density pc = 0.6po (in units of nuclear saturation density, po = 2.8 x 10^'* g cm~^), where 
nuclei merge into nuclear matter. The inner crust has ^ > Xc, where Xc = Rc/Rs > 0.9; 
for ^ > 0.9, the approximation -u = 1 — ^ is sufficiently accurate (Figure [T| left). The total 
momentum of inertia is 

/tot = ^^"!;^^ MP^ = 0.26MRI (1) 
while the inner crust has Mic = 4.9(1 — Xc^M and lie = 12.6(1 — Xc)^/i 



tot- 



The standard superfluid fraction, Q, is introduced to describe the protons and nuclei of 
the normal crust, Ic = (1 — Q)/tot, and the neutron superfluid component, Ig = Qhot, of the 
star. Although the neutron fraction varies with density, a typical average value is Q ~ 0.95 
( Zuo et aL]|2004 ). Regarding proton superconductivity in the core, here we choose the weak- 



drag limit. The model, however, gives reasonable results for the jump parameters also in the 
strong-drag limit; the decoupling of the core vorticity, pinned by fiux-tubes, reduces both 
ALgi and the moment of inertia responding at the glitch (cf. Paper I). 

The density dependence of /pin(p), the (mesoscopic) pinning force per unit length, is 
taken as in Figure [Tj^right). This is a reasonable first approximation for a parameter study 
in terms of the maximum value, f^', indeed, pinning goes to zero at Pc (no more nuclei) 
and at neutron drip pd = O.OOlSpo (no more neutrons), while it is expected to be maximum 
around densities where the pairing gap peaks. Moreover, we have performed a numerical 
simulation to evaluate /pin(p) in a bcc lattice, with random crystal orientations and proper 
vortex rigidity. We obtain profiles compatible with Figure [TJ with maximum values of order 



fm ~ 10^^ dyn cm~^ at densities pm ~ 0.2po. We also find, as already noted by Link (2009), 
that attractive and repulsive vortex-nucleus interactions are equivalent for pinning vortices 
to the lattice (Grill & Pizzochero, in preparation). 
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The geometry of the model is shown in Figure [2j The star spins around the ^-axis, and 
continuous vortex hues are assumed through the core. Following the results of |Ruderman &: 



Sutherland (1974), we can reduce the problem to axial symmetry by integrating the density- 



dependent quantities along the vortex lines; these quantities will then depend only on the 
cylindrical radius x = R/Rg, with R the distance from the rotational axis. We distinguish 
two cylindrical zones, separated by Xc = Rc/Rs- the 'crust' (x > Xc), with vortices lying 
entirely in the inner crust, and the 'core' {x < Xc), with vortices crossing the star core. In 
particular, we can integrate the pinning and Magnus forces to obtain an estimate of their 
total values on a vortex. If u{x) = Qs{x) — Q indicates the lag between the local superfiuid 
angular velocity and that of the rigid normal crust, the critical lag for depinning, Ci;cr(a^), is 
obtained by equating these two forces 



dz fpin[p{x,z)] =XUJcr{x)K / dzp{x,Z 



(2) 



where k 



h/rriN- 



In Figure |3| we show the resulting profile: UJcr{.x) presents a sharp peak in the 'crust', 
with maximum Wmax located very close to = 1 — Mm = 1 — Pml A; we will take pm = 0.2po- 
In most of the 'core', instead, ujcr{x) has a roughly uniform value Wmin ~ 10~^Ci;max (as x — )■ 
it diverges; similarly to the outer crust, however, this region can be neglected). We find (cf. 
Paper I) 



'^max ^criXrn) 



K Ad (7mag i^Xm ) 



(3) 



where 



2(1 -x) 



Vl — x^ — x^ In 



1 + v^ 



x^ 



X 



X 



vrx 



In 



1 + 



x" 



X 



x^ 



(4a) 
(4b) 



In the 'crust', this estimate of Ci;cr(x) should be reasonable, since pinning is continuous 

. In the 'core', instead, pinning 
is discontinuous: vortex lines are attached to the lattice only at their extremities, while most 
of their length lies in a pinning-free region (having selected the weak-drag limit). We can 
expect individual string-like excitations of the pinned vortices, which could detach them 
from the lattice well before cUmin is reached. Indeed, the collective rigidity of vortex bundles 
in coherent motion, which explains the axial symmetry predicted by the Taylor-Proudman 



along the vortices (Ruderman fc Sutherland|1974 Jones|1990 
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theorem ( Ruderman fc Sutherland||1974 ) , is actually not observed in laboratory experiments 



with superfiuid vortices attached to the rotating vessel only at their ends (Adams et al. 



1985). 



Although this issue requires and deserves further study, the crucial point is that vortices 
in the 'core' are pinned very weakly. On the other hand, drag forces due to magnetization 



correspond to very short relaxation times and thence very small steady-state lags (Alpar et 



al. 1984c). Moreover, Link (2009) has shown that vortex repinning is dynamically possible 
if the lag falls below a critical value (smaller than the critical lag for depinning). These 
considerations and the profile in Figure [3] naturally suggest the following scenario: as the star 
slows down, vortices in the 'core' are continuously depinned and then rapidly repinned; this 
dynamical creep allows a steady removal of the excess vorticity on short effective timescales 
Tc- Although the value of Tc is not relevant here, dynamical simulations of Vela glitches 



suggest Tc ~ 10° — 10^ s (Haskell et al. 2011), compatible with mutual friction dominated 



by vortex magnetization. On timescales At < Tc (e.g., during a glitch), the 'core' vorticity 
is only partially coupled to the normal crust (the rest being pinned or responding on longer 
timescales), and only a detailed study of the dynamics can provide a direct estimate of the 
coupled fraction. On longer timescales At 3> Tc, however, the dynamical creep ensures full 
effective coupling of the two components with (average) lag of order \fl\Tc; in steady-state, 
this scenario is then equivalent to the corotation model. 

The excess 'core' vorticity will be repinned in the 'crust', where pinning increases rapidly 
by orders of magnitude. At any time t after a glitch, the lag uj{t) = \floo\t defines a radial 
distance x(t) as in Figure [s] here floo indicates the steady-state (pre-glitch) angular acceler- 
ation. We now assume that the excess vorticity, corresponding to the entire region x < x{t) 
and to the lag w(t), is accumulated in a thin vortex sheet at x{t); as the star slows down 
and u{t) increases, the sheet is pushed outwards by the increasing Magnus force, and moves 
with x{t). When u(t) reaches the value cumax, the sheet is at the pinning peak, Xm, and 
the vorticity accumulated in the sheet is finally released simultaneously, causing the glitch. 
This picture reminds of a snowplow, pushing accumulated snow up an incline and eventually 
reaching its top edge. The interval between glitches is 

Atgi = (5) 

I ' 'oo I 

This scenario is compatible with post-glitch relaxation; indeed, after a glitch, the unpinned 
vortices in the 'crust' are under the same conditions as those considered in current post-glitch 
models. 



Although the 'snowplow' model is quite schematic, it contains a plausible mechanism 
for storing and releasing vorticity, as actually confirmed by parallel dynamical simulations 
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(Haskell et al. 2011). In particular, the model allows to calculate directly the angular mo- 

In Figure 111 we show the reduction of angular 



mentum Ly{x) of the vortex sheet at x 
momentum, = Ly{x)/L^{0), when uniformly distributed vorticity contained within x 

is accumulated in a sheet at x; at the peak, Xm, the reduction is of order 10~^. For com- 
parison, we also show the significantly different results for a uniform-density, cylindrical or 
spherical star; we see how spherical symmetry and realistic density profile are both crucial 
to obtain the correct order of magnitude of ^^(x). 

The angular momentum stored during Atg\ and released at the glitch, ALgi, can be 
calculated from the number of vortices removed from the interior and accumulated at Xm, 
namely AN^{xm) = Svri^^x^co'max/'t- We find (cf. Paper I) 



with an effective moment of inertia 



Iy(x) 



2(7r2 - 6) 



gv{x)QItot = TT g^{x)QMR^ 



(6) 



(7) 



where 



X 

~6 



2 r 



Vl-x^ {l + 2x^) -3x2 In 



1 + 



x^ 



X 



The glitch rise-time is very short, Tgi < 40 s (Dodson et al. 2002); we introduce a new 
parameter, Fgi, which globally describes the fraction of vorticity coupled to the normal crust 
on timescales of order Tgi (the steady-state coupled fraction, corresponding to long timescales 
and to pre-glitch conditions, is Y^o = 1). The value of Yg\ depends on the detailed short- 
time dynamics of the 'core' vorticity; in order to get an estimate of the observables, only 
this quantity is needed. From angular momentum conservation and variation of the crust 
equation of motion we find the glitch jump parameters 

ALgi 



/tot[l-Q(l-ygi)] 

Q(i - ^gi) 
i-Q(i-rgi)- 



(9a) 
(9b) 



3. Results and observations 

After fixing the basic stellar parameters M, and Q (more generally, M and an EoS), 
the model has two free parameters, fm and Ygi. It must predicts three observables: the 
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interval between glitches, and the jumps in angular velocity and acceleration during a glitch. 
In the case of Vela, the average observed values are Atgi ^ 3 years and Af2gi = 1.2 x 10~^ 



Hz (Lyne et al. 2000); we already mentioned that early observations gave Afigi/fioo ~ 10 



-2 



More recent data, however, indicate much larger values; in particular, the year 2000 glitch 



(Dodson et al. 2002) added to the already known short-, middle-, and long-time relaxation 
components (r^ ^ 10^ 10^ 10^ s, with AQi/Ct^ ^ 0.44,0.044,0.009 for i = 1,2,3), a fourth 
and very short one, with T4 = 1.2 ± 0.2 minutes and ACl4^/Cloo = 18 ± 6 (one sigma errors). 
In the 2004 glitch, however, such a component was observed only barely above noise and no 
firm conclusion could be drawn from the weak data ( Dodson et al.|[2007 ). Waiting for future 



observations, there is nonetheless evidence that right after a glitch Afigi/^oo is larger than 
unity. 

In order to test the model against observations, we consider a standard neutron star 
with M = 1AMq,Rs = 10 km, and Q = 0.95. If we take = 1.1 x 10^^ dyn cm'^ 
from equations [3]-[8] we find that Wmax = 0.01 Hz, and thence Atgi = 3.1 years and ALgi = 
9.5 X 10^^ erg s (also, AEgi = 6.7 x 10^^ erg). If we then take Ygi = 0.05, from equation 
9 we obtain Af2gi = 1.3 x 10~^ Hz and AQgi/Cloo = 9.3, in good general agreement with 
observations. In Paper I we analyze the parameter dependence of these results; we find that 
the model is quite robust under physically meaningful variations of all the basic parameters 

(M, Rs,Q,Pc,Pm,Pd)- 

In conclusion, assuming continuous vortices throughout the star, we find that maximum 
pinning forces of order fm ~ 10^^ dyn cm~^ can accumulate ~ 10^'^ vortices in the inner 
crust of a standard neutron star, and release them every ^ 3 years, transferring an angular 
momentum ALgi ^ 10'^^ erg s. This is one order of magnitude smaller than what inferred 
from the (microscopically inconsistent) assumption of disconnected vortices. Yet, it yields the 
observed glitch parameters, provided one assumes a small coupled fraction Ygi < 10%. The 
model is compatible with post-glitch recovery and with the presently known microphysics; the 
numerical results follow from implementing both spherical geometry and a realistic density 
profile, and they are robust. 

This work was supported by CompStar, a Research Networking Programme of the 



European Science Foundation (http://www.compstar-esf.org/). 
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Fig. 1. — (left) Density profile for the n=l polytrope in dimensionless variables; the dotted 
line is the low-density approximation m = 1 — ^. In the insert, the low-density region is 
magnified, (right) Density dependence of the pinning force. 
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Fig. 2. — Geometry of the model (not in scale): the star rotates around the 2;-axis, while the 
X-axis represents the dimensionless radius in cylindrical geometry. Two vortices are shown 
at xi and X2, respectively in the 'core' and the 'crust', the two cylindrical regions separated 
by Xc = Rc/ Rs ■ 
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Fig. 3. — Critical lag for depinning, uJct{x), as a function of the dimcnsionless cylindrical 
radius x. At any given time, the lag u!{t) determines the position x{t). The lags have been 
normalized to cUmax- 
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Fig. 4. — Reduction of angular momentum, £v{x) = L^{x) / L,^{{)), when uniformly distributed 
vorticity contained within x is accumulated into a vortex sheet at x. Three scenarios are con- 
sidered: uniform-density cylindrical star (dotted), uniform-density spherical star (dashed), 
n=l polytrope (solid). In the insert, the 'crust' region is magnified. 



